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Abstract 

Non-equilibrium quantum field theory studies time dependence of processes which are 
not available for the S-matrix description. One of the new methods of investigation in non- 
equilibrium quantum theory is the stochastic limit method. This method is an extension of 
the works by Bogoliubov, van Hove and Prigogine and it permits to study not only the system 
but also the reservoir degrees of freedom. We consider the stochastic limit of translation 
invariant Hamiltonians in quantum field theory and show that the master field satisfies a 
new type of commutation relations, the so called entangled (or interacting) commutation 
relations. These relations extend the interacting Fock relations established earlier in non- 
relativistic QED and the free (or Boltzmann) commutation relations which have been found 
in the large N limit of QCD . As an application of the stochastic limit method we consider the 
photon splitting cascades in magnetic field and show that photons in cascades form entangled 
states ("triphons") and they obey not Bose but a new type of statistics corresponding to the 
entangled commutation relations. 
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1 Introduction 



The basic object to study in quantum field theory is the S-matrix introduced by Heisenberg. 
Bogoliubov and Shirkov have developed the S-matrix formalism which includes all the quantities 
considered in quantum field theory |T[. The physical idea behind the S'-matrix approach is that in 
the scattering processes there exists a characteristic time scale such that in a time regime larger 
then this time scale one can neglect interaction and particles evolve according to the free dynamics. 

The situation in statistical physics is different because here one has not just one but several 
relevant time scales and as a result we don't have here a universal method comparable with the 
S'-matrix approach in quantum field theory. One can say that the role of S-matix approach in 
non-equilibrium statistical physics is played by various master and kinetic equations. It was the 
fundamental Bogoliubov idea about the existence of two time scales which lead to the modern 
progress in the microscopic derivation of kinetic equations || . 

Methods of quantum field theory, in particularly Green functions, are widely used in equilib- 
rium and non-equilibrium statistical physics [§, J|, ||. It is well known that the so-called double- 
time thermodynamic Green functions which were used by Bogoliubov and Tyablikov H have had 
great success especially when applied to magnetic problems. From the other side there are not 
so many investigations in quantum field theory where results and methods of n on- equilibrium 
statistical physics are exploited. N.N. Bogoliubov has made great contributions to both of these 
sciences. As far as we know his works in non-equilibrium statistical physics and in quantum field 
theory were performed completely separately. Probably the reason was that important problems 
in quantum field theory in that time were related with the scattering processes described by the 
S-matrix and they were really very different from typical problems in non-equilibrium statistical 
physics such as the derivation of kinetic equations. 

We would like to point out that there are important problems in quantum field theory where 
the standard S-matrix description is not very convenient or even not applicable. These include not 
only investigation of bound states and spectral problems (see 0) but also processes with unstable 
particles [0, |8| (in fact almost all particles are unstable), atom-photon interactions H, elementary 
particles in "semidressed states" with non-equilibrium proper fields [JTO]] , electroweak baryogenesis 



and phase transitions in the early Universe and in high-energy collisions ||11||, quantum optics 



12f etc. In the consideration of such processes we are interested in the time regime smaller than 



the "infinite" time when the S-matrix description becomes applicable. One can say that the 
consideration of such processes belongs to non- equilibrium quantum field theory. One believes 
that the method of S-matrix in quantum field theory is analogous to the Gibbs distribution 
in equilibrium statistical physics and that there exists a general method (the stochastic limit, 
see below) in non-equilibrium quantum field theory which provides a description of quantum 
phenomena depending on time. One of the first works on the systematic application of methods 
of non-equilibrium statistical physics in quantum field theory is the work of Prigogine [13|] in which 
kinetic equations for the Lee model have been derived. 

A general method in non-equilibrium quantum field theory is the method of stochastic limit 
rjL ITH|. The idea of this method is the systematic application of the A 2 t-limit and quantum 



stochastic differential equations. One considers the evolution operator U(t) of quantum system 
for small coupling constant A and large time t. Intuitively, the weak coupling - long time limit 
means that we are looking at times in which the particle has already weakly interacted many times 
with the field (long time cumulative effects). The net average effect of these interactions amounts 
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to a loss of memory (Markovian approximation). Therefore in this limit we can expect to be able 
to approximate the microscopic time evolution, which contains complicated memory effects, with 
a simpler markovian evolution. 

When we say that the S"-matrix method is not sufficient in non-equilibrium quantum field 
theory we mean that the standard dynamical definition of the S"-matrix in real time is given 
for example in terms of wave operators (including dressing [T7|, [TJfl) or LSZ-formalism. This 



definition can not be applied immediately to the processes with unstable particles. The flexible 
Bogoliubov- Shirkov approach to S-matrix in principle can be applied to the description of 
unstable particles. There exists a phenomenological approach to S'-matrix which is not based on 
a Hamiltonian formalism. In this approach unstable particles are described by the Breit-Wigner 
complex poles of the scattering amplitudes |19[]. The dynamical justification of this phenomenolog- 



ical approach is given in the Weisskopf-Wigner resolvent method, for a discussion of the resolvent 
method see for example |8], The resolvent method is usually used for the investigation of 
the degrees of freedom of the system interacting with reservoir. The stochastic limit method is 
conveniently used for the consideration of degrees of freedom not only of the system but also of 
the reservoir. 

The first rigorous result about the interaction of a system with a reservoir where the role of 
a scaling limit involving \ 2 t begun to emerge, is due to Bogoliubov [fZI]] (see the next section). 



Friedrichs, in the context of the now well known Friedrichs model p2| , was lead to consider the 
scaling limit 

A — > 0, t — > oo, \ 2 t = constant (1) 

by second order perturbation theory. 

In the mid 1950's van Hove used this scaling as a device to extract the effects of a small 
perturbation of the global Hamiltonian of a composite system, on the reduced evolution of a 
subsystem and to derive a Pauli-type master equation describing the irreversible time evolution 
of the observables of the sub-system. Therefore, in the quantum theory of open systems, the limit 
([]]) is known as the van Hove or the \ 2 t limit. 

Using the perturbative development of the dynamics in powers of A, van Hove argued that the 
terms of order 2n should behave as A 2n t n (in contrast with the rough estimate X 2n t 2n ). Furthermore 
he made plausible that almost all the terms of order 2n should behave as X 2n t n ~ £ , for some e > 
and therefore vanish in the limit ([!]) , while the remaining terms should sum so to give a transport 
equation. 

In the stochastic limit we study the behaviour of the system in the large time and small coupling 
constant regime. If G(x\, ...x n ) is the Green function then we want to investigate the asymptotic 
behaviour of the expression G{xi, x\/\ 2 , x n , x°/A 2 ) when A — > 0. This can be performed by 
using the anisotropic renormalization group method [I, [23], [23 1. 



Notice that the A 2n t n -behaviour of the terms of order 2n of the iterated series is exactly of the 
same order of magnitude of the behaviour of the moments of order 2n of a Brownian motion with 
time parameter rescaled as X 2 t. Therefore a posteriori we can interpret van Hove's perturbative 
result as a fist indication that the limit ([!]) should evidentiate a kind of quantum Brownian motion 
or quantum white noise underlying the dynamics of the quantum system. To describe and to solve 
the dynamical equations after the stochastic limit one has to derive the stochastic limit for the 
collective operators, the so called master field. For simple models the master field is the quantum 
white noise whose creation and annihilation operators satisfy the relations 

[b(t)X(t')]=S(t-t') (2) 
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In the stochastic limit one gets for the evolution operator the equation 

dU(t) 



dt 



(F + (t)b(t)+b + (t)F(t))U(t) 



Here the white noise operators b(t),b + (t) are singular operator functions of t and F(t) is a reg- 
ular operator function of t. This singular equation represents the Hamiltonian form of quantum 
stochastic differential equations. This equation can be explicitly solved for many models. The 
crucial role for the possibility of the explicit solution plays the commutation relation ([!]). 

For more complex models the master fields are more complex and one gets the entangled 
or interacting commutation relations. The program of investigation of models of quantum field 
theory in the stochastic limit consists from two parts. First we have to find the commutation 
relations for the master field and then study the singular differential equation for the evolution 
operator. In this paper we shall discuss only the first part of the program. For the consideration 
of the evolution operator see for example |T5L E31 



In recent years various modifications and deformations of the algebra of canonical commutation 
relations have been discussed, see for example [p7| [28], and references therein. In particular, 
in the large N limit for SU(N) invariant gauge theory (as well as for NxN matrix models), the 
following relations 

b{k)b + {k') = 5{k - k') (3) 



appear naturally [E]J. Here k, k' are momentum variables. The algebra generated by the operators 
b(k),b + (k') satisfying ([3D is called the free (or Boltzmann) algebra. 

In this paper we prove that the stochastic limit of interacting fields, under the only constraint 
of momentum conservation, leads to a new algebra of commutation relations. We get the relations 
of the form 

Biplh, k m )B + (p'\k[, k'J = n(p)5(E(p, h, k m ))5(p - p')8{k x - k[)...5(k n - k' n ) (4) 

where B (p\ki, k m ) is the master field, obtained as the stochastic limit of a translation invariant 
interaction Hamiltonian, n(p) is the operator density of particles, and E(p, k±, k m ) is the energy 
associated with the interaction vertex. These relations are called the entangled (or interacting) 
commutation relations. Notice that the equality (Q) extends the Boltzmann algebra ([3]) because 
the right hand side is an operator (in the particle space) rather than a scalar: in this sense one 
speaks of Hilbert module |3(J (rather than Hilbert space) commutation relations. B(p\ki, ...,k m ) 



and the density n(p) satisfy 

[n(p'), B(p\k u k m )\ = (6(p' -p)- 6{p' -P + J2 h))B{p\k u k m ) (5) 

As one of physical applications of the above ideas we will argue that photon splitting cascades 
in the magnetic field create entangled states and that photons in cascades obey not Bose but a 
new type of statistics - infinite or quantum Boltzmann statistics. Therefore, this statistics has a 
physical meaning since it describes photons in cascades and more generally the dominant diagrams 
in the long time/weak coupling limit in quantum field theory. The states in cascades are formed 
from triples of entangled photons and may be called triphons. They belong to an interacting Fock 
space |3(], |3l[] . Interacting Cuntz algebra has been considered in p2| . 



The standard definition of the stochastic limit is given as the limit of the Wightman correlation 
function. For some models the limit of these correlation functions is equal to zero. We show that 
the stochastic limit for the Green functions is non-trivial even in these cases. 
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The paper is organized as follows. In Sect. 2 we remind the Bogoliubov theorem and compare 
it with the stochastic limit. In Sect. 3 we discuss the stochastic limit. Sections 4 and 5 are devoted 
to the derivation of the entangled commutation relations. In Sect. 6 we discuss the stochastic 
limit for the Green functions and in Sect. 7 the universality classes of the stochastic limit and 
decay processes are considered. Finally Sect. 8 is devoted to applications of the stochastic limit 
to photon splitting cascades. 



2 Bogoliubov's Theorem 

The first rigorous result, about the interaction of a system with a reservoir, where the role of a 
scaling limit involving A 2 t begun to emerge, is due to Bogoliubov |21| who considered a classical 
system with Hamiltonian 

H = H S + H R + H mt (6) 

where 

j ^ N N 

H s = -(p 2 + u; 2 q 2 ), H R = -J2{pl + u 2 n q 2 n ), H int = Ag £ (a n q n ) (7) 

A A 71=1 71=1 

He assumed that J2<r<u> n a 2 /w 2 ~~ > XT' J{y)dv < oo as N — ► oo with J(v) being a non-negative 
continuous function. Supposing that, at t — the system is in the state q = q ,p = p and that 
the state of the reservoir is a random variable with the distribution p R = exp{—H R /T} he proved 
that: 

• As iV — > oo the limit distribution p s = ps(t,p, q) of the random variables p = p t ,q = q t , of 
the system, exists at any time t > 0. 

• If Hs is given by (2), define the function p° s = p (t, H s ) = pa(t,p, q), by: 



p°s(t,H S ) = 4?r2T(1 ! e - 2aA2f) J /; #e XP {- — — r(1 _7-^f — ^} 

where a = jJ{oj) and E is (pg + <Zo)/2 where (po, Qo) is the initial state. Then for small A 
the function ps can be approximated by the function p° s , in the sense that, for any positive 
a, (3 as A — > 0, one has, uniformly in the interval < t < 



2tt 



H s + E e- 2aX ' 2t - 2^Hs~E^e~ 2aXH cos, 



l r t+At >- 

^ J (Ps - P° s )dt - 



for any subsequence {At a} such that A 2 At a — » 0, At a — > c>o 

He gave a rigorous proof of this result using a Volterra integro-differential equation. 

Bogoliubov's condition A 2 AtA — > is different from the one used in the stochastic limit (A 2 t — > 
constant), notice however that p s (t) depends only on A 2 t and also that for t — ► oo one gets the 
Gibbs distribution, 



-e 



-H S /T 
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3 The Stochastic Limit 



The stochastic limit is now widely used in the consideration of the long time/weak coupling 
behaviour of quantum dynamical systems with dissipation, see for example ||15|| . 
Let be given a quantum system described by the Hamiltonian 

H = H + XV 

where A is the coupling constant. The starting point of the stochastic limit is the equation for the 
evolution operator in interaction picture 

where 

V(t) = e itH *Ve~ itHa 

The main idea is that there exist a new quantum field (master field) and a new evolution operator 
U(t) (they both live on a space different from the original one) which approximates the old one 

U w (t) «W(A 2 t) 

and the approximation is meant in the sense of appropriately chosen matrix elements. The above 



approximation suggests a natural interpretation of the van Hove rescaling [16 A — ► 0, t — > oo 
so that X 2 t ~ constant = r (new time scale): it means that we measure time in units of 1/A 2 
where A measures the strength of the self-interaction. By putting r = 1 we see that the van Hove 
rescaling is equivalent to the time rescaling t — > t/X 2 , and therefore the limit A — ► will capture 
the dominating contributions to the dynamics in the new time scale (the error can be estimated to 
be of order A 2 ). It is remarkable that in this limit the dominating contributions can be explicitly 
resummed giving rise to a new unitary operator. 

A simple change of variables shows that the time rescaling t — * t/X 2 is equivalent to the 
following rescaling of the Schrodinger equation for the evolution operator: 

= -i{v W xmm (8) 

The unitary operator U(t) is then obtained by taking the limit A — > 0: 

U(t) = \imU w (t/X 2 ) (9) 

A — »0 



and the corresponding limit equation is 

dUjt) 
dt 

where 



-iV(t)U{t) 



-v(± 

a-oA VA 2 

For a number important models the interaction Hamiltonian has the form 



V(t) = hm T V{- 



V = A + JC 
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where A is a monomial in the creation and annihilation operators. The master field is given by 
the asymptotic behaviour of the collective operator 

and its Hermitian conjugate. Here 

A(t) = e UHo Ae- itHo 

The stochastic limit is meant in the sense of the convergence of correlation functions 
< A ei (t 1 /X)...A en (t n /X) >. Here e, = ±. 



4 Entangled Commutation Relations 

In recent years there has been an interest in various modifications and deformations of the alge- 
bra of canonical commutation relations. Especially the simple relations (H) were discussed. An 
extension of the algebra (|3|) has been given in [J3TJ] as the algebra describing the interacting Fock 
space |3(| obtained in the stochastic limit for non-relativistic QED. 



In this paper we prove that the stochastic limit of interacting fields, under the only constraint 
of momentum conservation leads to a generalization of the algebra (§). We obtain that the new 
algebra has as its generators the master field B(p, k) depending on two momenta p and k and the 
operator density of particles n(k) which satisfy the relations 

B(p, k)B + (p', k') = n{p)5{E{p, k))5(p - p')5{k - k'), (10) 
[n{p'), B(p, k)} = W - p)) - 6(p' -p + k))Bip, k), (11) 

[n(p),n(p')]=0 (12) 

Here E{p, k) is the energy associated to the interaction vertex. 

We call the relations (|T0"|)-(|r2"|) entangled commutation relations because, on one hand they 
allow to calculate correlations of any order among the field and, on the other hand they show that 
the master fields are not kinematically independent. 

In the construction of these operators we use the Van Hove time rescaling t — > t/\ 2 , A — > 
where A is the coupling constant. Q 

One can get a generalization of the algebra (jT0"|)-(|T2"|) for the multiparticle master field B(p\k\, k n ) 

In [ 31 1 a generalization of the algebra (|j) of the following form has been obtained: 



Bik)B + ik') = 5ik - k')5iuik) + P ■ k) , [P,Bik)] = kBik) (13) 
with P being the operator of momenta of particles. This algebra is not realized in the usual Fock 



space but in the interacting Fock space. One can get the algebra (p~3|) from the algebra ([10] 
if we set 

P — I pnip)dp, Bik) = / Bip,k)dp 
lr To avoid a discussion of renormalization procedure we assume that there is an ultra-violet cut-off 
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These relations show that, contrarily to what happens before the limit, the observables of the 
particle do not commute with the master field. In other words: before the stochastic limit the 
particle and the field are kinematically independent but nonlinearly related by the dynamics; after 
the stochastic limit the dynamics is simplified but particle and field are no longer kinematically 



independent: this is what we call entanglement. These new features imply that the algebra (10) 
flT2] ) is not realized in the usual Fock space but in the interacting Fock space [3C|. 
Let us consider the Hamiltonian 

H X = H + XV 

where the free Hamiltonian 

H = f e(p)c + (p)c(p)d 3 p+ f uj(k)a + (k)a(k)d 3 k, 



{c(p), c + (p')} = 5(p - p'), [a(k), a + (k>)] = 5(k - k') 
and the interaction Hamiltonian: 

V = [ d 3 kd 3 pg(k,p)(c + (p)c(p - k)a(k) + h.c.) (14) 



Here g(k,p) is a test function and e(p) and u(k) are one-particle dispersion laws, for example 
e (p) =p 2 /2, oo(k) = \k\. 

The rescaled collective fields in this case have the form 



A x (p, k,t) = \ e^c + {p)a{k)c{p - k)e~^ = ~ c + {p)a{k)c{p - k) e itE{p ' k)/x2 (15) 
A A 

A+{p, k, t) = - e^c + {p - k)a + {k)c{p)e~^ = - c + {p- k)a + {k)c{p)e- ltE(p ' k)/x (16) 

where 

E(p, k) = e(p) — uj(k) — e(p — k) (17) 
is the corresponding energy. One has the following main theorem. 

Theorem 1. The stochastic limit 

lim A x (p, k,t) = B~(p,k,t), lim At (p,k,t) = B + (p,k,t) 

A- — >ij A — 

exists in the sense of the convergence of the matrix elements ( q = ±) 

lim < 0\c{q)A x 1 {p 1 ,k 1 ,t 1 )...A x n {p n ,k n ,t n )c + {q')\0 >= 

A — ^0 

(^ ,c(q)B^( Pl ,k 1 ,t 1 )...B £n (p n ,k n ,t n )c + (q')^ ) (18) 

as distributions and the limiting operators B~ = B and B + satisfy the entangled commutation 
relations 

B(p, k, t)B + (p', k', t') = 2n5(t - t')5(p - p')5(k - k') ■ 5(E(p, k))n{p) (19) 
[nip'), B*(p, k, t)} = (±)(5(p' -p)- dip' -p + k))B*(p, k, t) (20) 



S 



[n(p), nip')} =0 (21) 

Here \l/o is the vacuum in the new Hilbert space, B(p,k,t)Y\iC + (qi)^o = 0. We use the same 
notations for the creation and annihilation operators of c- particles in the original and in the new 
Hilbert spaces. n(p) is the operator density of the c-particles, n(p) = c + (p)c(p). 
If we set 

B(p,k,t) = b(t)®B(p, k) 

where 

b{t)b + {t r ) =2ir5(t-f) 

then we get the relations ([10) 

B(p, k)B + (p', k') = n(p)5(E(p, k))5(p - p')5{k - k') 

The theorem will be proved in the next section. Note that to get non-zero in the RHS of (|19|) 
we have to chose suitable dispersion relations, so that there are non-trivial solutions of equation 
E(p, k) = 0. We will discuss this point in Sect. 6. 

5 Proof of the Theorem 

Let us consider the matrix elements 

(0|c(g)n^fe,M*)c + (g')|0> (22) 

i=l 

To evaluate ( [22] ) we apply Wick's theorem and consider only the corresponding connected 
diagrams. Each vertex contains 3 lines which are characterized by two momenta (ki,pi). We find 
the momentum corresponding to the 3rd line using the momentum conservation. If a diagram 
contains L loops then only L + A + B — 1 momenta are independent (L momenta for loop variables 
and A + B — 1 momenta for the exterior lines). 

For every vertex there is the corresponding energy exponent. These energies depend on the 
momenta of the lines that enter to the given vertex, 

Ef = Et(k llPl ) 

via dispersion laws, 

Ef = [e(pi) - e( Pl ± h) ± u{ki)) (23) 
In the proof of the theorem we will use the following lemma. 

Lemma 1. One has the following relation in the sense of distributions 

\iin^e itE{p ' k)/x2 = 2rr5(t)5(E(p, k)) 



The proof of the lemma is standard, see [3~3| 
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i 



Figure 1: Diagram with a pair of conjugated vertices 



Let us consider a connected diagram corresponding to the matrix element fl22|). The proof 
of the theorem consists of three parts. First we prove that only diagrams consisting of pairs of 
conjugated vertices don't vanish in the limit A — > 0. Next we prove that such diagrams are in 
fact non-crossing or half-planar diagrams. And finally we show that the non-crossing diagrams 
are described by the entangled commutation relations. 

Generally, the sets of momenta corresponding to different vertices are different. However, it 
may happens that the same set of momenta corresponds to two different vertices. More precisely, 
momenta which come in the first vertex come out from the second one and viceversa. 

Definition 1. We say that two incident vertices of a given connected diagram are conjugated 
if the momenta coming in the first vertex come out from the second vertex, i.e. the vertices have 
the same momenta but with the opposite orientation. 

If the i— vertex has a conjugated vertex then we denote the latter by i. A typical example of 
diagrams containing at least one pair of conjugater vertices is a diagram with a mass insertion 
such that this insertion contains a line that does not cross others lines of the diagram (see Fig.l). 

Let us prove the main 

Lemma 2. // a connected diagram doesn't consist only from pairs of conjugated vertices then 
it vanishes in the limit A — ► (in the sense of distributions). 

Proof. To a given diagram, representing a matrix element (^) being integrated over t±, . . . , t v 
with test functions, corresponds the expression that schematically can be written as 



here by t we mean t\, . . . ,t v , by p we mean pi, . . . ,pa, p[, . . . ,p' B and q denotes the set of inde- 
pendent momenta associated with the diagram under consideration. Ei are given by (^) and 
4>(t,p,q) is a test function. 

To evaluate the asymptotic behaviour of this expression when A — > we will make the change 
of variables corresponding to the conjugated vertices. Notice that for the diagram doesn't vanish, 
the number of vertices v should be even, v = 2n. Suppose that there are no pairs of conjugated 
vertices which are denoted ...,i no ,i no }. Let us divide the set of all vertices {1, ...,2n} into 




1 



A+B-l L 




(24) 



a 1=1 i=l 
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two disjoint subsets i 2 , . . . , i no ,i no+ i, i n } and {h, i 2 , . . . , i no , in+no+x, ■ ■ ■ hn} such that in 
every subset there are no conjugated vertices. We denote the corresponding set of time variables 

{th, -tin} by t{1) and the set { t i 1 »"ii„/i„ + „ 0+1 ,"4J b y t{2) - 

Now we perform the following change of variables 

(t«t( 2) )^M 2) ) (25) 

t (D = t (2) + A 2 T; 

or more precisely 

(t«i( 2) ) = (ti,...,t 2n ) -> (r,t {2) ) = (n,...,r n ; tg.,j = 1 . . . ,n ; £ in+r , r = n + 1, . . . ,n) (26) 

fc. = ^ + \ 2 Tj , 1 < j < n (27) 

^ = t in+ . + X 2 Tj , n <j <n (28) 
After this the integral takes the following form 

*E; =1 ^ e *E^+V^ /A2 . *E; = „ 0+1 (^+^ n+i )W* 2 



0(t( 2 )+A 2 r,t( 2 ), P ,g)-n^n^ j n ( 29 ) 

j=l j=l j=n +l 

By definition of conjugated vertices E{ . + = and we left with 



/" iV" iV" . +Ej , .)ti , ■ / A 2 



0(t( 2 )+A 2 r,t( 2 ), P ,g)-n^n^ n ^n^-n^ ( 3 °) 

j=l j=l j=n +l 

Here t^ 2 ^ + A 2 r schematically represents the dependence of the half of the t-variables on A. When 
A — > we can neglect the dependence of (ft on A and the integration over r gives the product of 

/n*(^) ci ^ +i( ^ + ^ iWA V(*,*,p,ff) n ^n*n d 9 ( 3i ) 

3=1 i=«o+l 

Note that the second exponent in the expression ([!(]) vanishes since the energies in the conju- 
gated vertices are equal. 

Suppose that in the diagram there are non-conjugated vertices, i.e. n ^ n Q . When A — > the 
expression (^1|) goes to zero since, according to our assumption the set of momenta in vertices 
ij and i n +j, no < j < n do not coincide and therefore the functions E^ + E in+] as functions of 
momenta don't vanish and therefore according the Riemann-Lebesgue lemma we get zero in the 
limit A -> 0. 

In the case when n = n the exponent in fl3T|) vanishes and generally we get the non-zero 
answer: 

/ (fi Wj <P({h}A%hP,<l)I[dpl[dq (32) 
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The lemma is proved. 

Lemma 3. If a connected diagram with two external lines consists only from pairs of conjugated 
vertices then it is half-planar, i.e. it can be drown in the half-plane without self intersections. 

We will not present here the simple proof of this lemma. 

Lemma 4. The limiting expression is equal to the RHS of the relation 

The proof is obtained by using the algebra ([HJ), (|2"T|). 

Now the theorem follows from the above four lemma. 

Remark 1. We have considered only matrix elements with two external particles and as a 
result the matrix element doesn't vanish only if v is even, v = In. In the general case of matrix 
elements with an arbitrary number of external lines the number of vertices v can be odd. For 
the case of odd v, v = 2n — 1, we once again select n vertices and make the following change of 
variables 

(t 1 ,...,t 2n -i) -> (r 1 ,...,T n ;ti.,j = l,...,n ;t in+r ,r = 1, . . . , n - 1 - n ) (33) 
with the Ti as before (see (f2~7|)). 

Uj = % + A 2 r 3 - , 1 < j < n 
U 3 = U n+j + \ 2 Tj n + l<j<n-l 

The difference with the case of even v is that we get an extra factor A in front of the integral 

n no 

<f)(t + Ar, t, p, q) Y[ dpdq J[ dtj J[ dt in+j (34) 

j=l j=n 

Here we use the same schematical notations as in (5). When A — > we neglect the A-dependence 
of <p and we get a product of 5-functions. The second exponent goes out. The third exponent 
disappears in the case then no = n — 1, but since we have an extra factor A the expression (8) 
always goes to zero as A — > 0. 

Remark 2. The theorem admits a generalization to the case of an arbitrary number of external 
particles 

< 0\]Jc(q l )Al 1 (p 1 ,k 1 ,t 1 )...Al n (p n ,k n ,t n )Y[c + (q' J )\0 > 

i j 

Remark 3. For the composite operators 

A x (t,p, k u ...,k m ) = ^c + {p)c{p-k 1 )a^{k 1 )...a^^^ 
in the stochastic limit one gets the entangled commutation relations (£§). 



It). 
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6 The stochastic limit for Green functions 

The master field in the standard formulation of the stochastic limit is defined by the limit at 
A — > of the Wightman correlation functions 

<0\A x 1) (t l )...A x n \t n )\0> (35) 

This limit defines the Hilbert space in which the master field lives. For some models the limit of 
these correlation functions is trivial (equal to zero). However this doesn't mean that the stochastic 
limit of such the model is trivial because we can consider another natural formulation of the 
stochastic limit which is defined by the convergence of the chronologically ordered correlation 
functions (Green functions) 

< 0|T(4%)...4%))|0 > (36) 
For the evolution operator one has to consider 

<0\T(4\t 1 )...A?(t n )U(t))\0> (37) 
Here the T-product is defined as 

T(A$\t 1 )...A™(t n )) = At\t h )...4 n \K) 

if t h > ... > t in . 

The limit of the Green functions can be nontrivial even if the limit of the Wightman functions 
vanishes. In the simplest case in the first formulation we use the relation 

lim^e ite / A2 = 2nS(t)5(x) (38) 

and in the second 

H^^'-^Th (39) 

Here 6{t) = 1, t > 0; 6{t) = 0, t < 0. In particular if 

A x (t,p) = je u ^ x2 a(p) 
where [a(p), a + (p')] = 5{p — p'), p,p' G R 3 then the stochastic limit of the Wightman functions is 
< 0\A x (t,p)A+(t',p')\0 >= J-e^-'MpV^fc-p') -> 2n5(t - t')5(uj(p))5(p - p') 

A 

and the stochastic limit of the Green functions is 

< 0\T(A x (t,p)At(t',p'))\0 >= - fle^Mrt/^fr - p') - iS(t - t') * + J (p - p>) 

It is important to notice that in the last formula one gets a nontrivial limit even if 5(u(p)) = 
when in the former formulation the limit is trivial. 
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7 Stochastic Limit and Decay 

In the formulation of Theorem 1 we assumed the special form of the interaction Hamiltonian 
(|14T). In local quantum field theory the typical Hamiltonian is more complicated than flTil) . In 
particular, for the Yukawa interaction of fields ip and <fi the Hamiltonian has the form 

H X = H + XV 

where Hq is the sum of the free Hamiltonians for the fermionic field ip and for the scalar field <fi 
with relativistic dispersion laws, 



u a (k) = y/ml + k 2 , u b (k) = sjm 2 b + k 2 (40) 

and 

V = J (fxg^cp = 

d 3 kd 3 pg(k, p)(c + (p)c(p - k)a(k) + c + (p)c + (k - p)a(k) + c + (p)a + (k)c + (-p - k) + h.c.) (41) 



Here g{k,p) is a test function. This expression is not a well defined operator in the Fock space 
but it defines a bilinear form. We have the following collective operators 

1 it 

— c + (p)c(p — k)a(k) exp — (e(p) — e(p — k) — uo(k)) (42) 
A X 2 

1 it 

- c + (p)c + (k)a(p + k) exp — (e(p) + s(k) - u(p + k)) (43) 
A X A 

1 it 

- c + (p)c + (k)a + (-p - k) exp — (e(p) + e(k) + u(p + k)) (44) 
A A 



In the limit A — » the operator (f44|) vanishes because in the correlation functions one gets 

5(u c (p) + uj a (k) + to c (p + k)) = 

due to the positivity of energy. 

The limit of the operator (fE^) is also zero since for relativistic dispersion laws (|40|) 

S(u c (p) + u a (k) - uj c (p + k)) = 

Only operator (|43|) has a chance to be non-zero. In this case we have to find non-trivial solutions 
of 

LOcip) + LO c (k - p) = LO a (k) (45) 

There are solutions if 

m 2 > 2m 2 



i.e. if one has the decay. Therefore we obtain that the relativistic interaction Hamiltonian (|4lD is 
within the same stochastic universality class as the interaction Hamiltonian ( |T4D 

V = [ d 3 kd 3 pg(k,p)(c + (p)c + (k - p)a(k) + h.c.) (46) 
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In quantum electrodynamics (QED) the interaction has the form 



V = J tfxilrf^Ap 

If one neglects the spinor and polarization indices then the free Hamiltonian has the form ( |40|) 
and the interaction Hamiltonian has the form (JO]). We see that the stochastic limit of QED 
Hamiltonian is reduced to the form of the above discussed Hamiltonian and as a result we obtain 
the trivial limit, since there is no decay in the standard QED. However, if we consider QED in 
external field we can get a non-zero result due to the change of dispersion law. One of such 
examples we will consider in the next section. 



8 Photon splitting cascades and a new statistics 

In the previous section we have seen that the stochastic limit of Wightman functions in relativistic 
QED is trivial and that it can be non-trivial in the presence of external fields. It is known that 
there is the splitting of a photon into two in an external magnetic field. This splitting is one 
of the most interesting manifestations of the non-linearity of Maxwell's equations with radiative 
corrections in an external magnetic field. In a constant uniform field, this process occurs with 
conservation of energy and momentum. The process was considered by Adler et al. in the early '70s 
by using the Heisenberg-Euler effective Lagrangian |34], |35|, |36| . Photon splitting was considered 
as a possible mechanism for the production of linearly polarized gamma-rays in a pulsar field. 
Recently the splitting of photons has found astrophysical applications in the study of annihilation 
line suppression in gamma ray pulsars and spectral formation of gamma ray bursts from neutron 
stars |37|, ill- Photon splitting cascades have also been used in models of soft gamma-ray repeaters, 
where they soften the photon spectrum |$9], [EJ. The process of photon splitting is potentially 
important in applications to a possible explanation of the origin of high energy cosmic rays from 
Active Galactic Nuclei [^TJ. A recalculation of the amplitude for photon splitting in a strong 



magnetic field has been performed recently in [[42 



We will start from a discussion of the theory of photon splitting cascades and show the emer- 
gence of infinite statistics in this theory and then discuss its connection with the stochastic limit of 
quantum field theory. In the decay of a photon with momentum k into photons with momentum 
ki and k 2 , we have conservation of momentum and energy k = k\ + k 2 , u>(k) = oj\{ki) + ^(A^). 
For photons in vacuum, in the absence of external fields, uj = uoi = uj 2 = k and although these 
two equations have a solution the decay is forbidden by the invariance under charge conjugation 
(Furry's theorem). 

In a constant uniform magnetic field Bq there are only two decay processes kinematically 
allowed, 7y — > 7_l+7_i_ and 711 — > 7||+7_i_ p5| . Here the subscripts _L and || will denote polarizations 



of the photon with respect to the vector B . More precisely, in presence of a magnetic field one 
has a distinctive plane, namely the kB plane. One takes the linear polarization of the magnetic 
field of the photon parallel and orthogonal to this plane as the two independent polarizations of 
the photon, || and _L, respectively. 

The vacuum in the presence of the field Bq acquires an index of refraction n, and the photon 
dispersion relation is modified from k/uo = 1 to k/uo = n. The indices of refraction nn± can be 
calculated from the Heisenberg-Euler effective lagrangian. Adler showed that for subcritical fields 
in the limit of weak vacuum dispersion only the splitting mode || + _L operates below pair 
production threshold. For weak dispersion n± = l+^j ft and n\\ = P, where (3 = Bq sin 2 9 
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Figure 2: Cascade: 7y(p) -> 7_l(& 1 ) + 7y(p - fci) -> 71 (&i) + 7±(fe) + 7||(p - &i - fc 2 ) 



and is the angle between /c and B . It is mentioned by Harding et al. that in magnetar models 



of soft gamma repeaters , where supercritical fields are employed, moderate vacuum dispersion 
arises. In such a regime, it is not clear whether Adler's selection rules still endure since in his 
analysis higher order contributions to the vacuum polarization are omitted. In [|38| photon cascades 
are considered for the case where all three photon splitting modes allowed by CP invariance are 
operating. Baier et al. have found that there is only one allowed transition (|| — + _L) for 
any magnetic field. They suggested that a photon cascade could develop only if magnetic field 
changes its direction. It seems that the question on the validity of Adler's rule for a non weak 
vacuum dispersion deserves a further study. In this work we consider photon cascades when both 
kinematically allowed modes (|| + _L and || — > ||+ _L) operate. 



The interaction operator for the decay || + _L is known to be [|36] 

V x {t) = X 1 J (B E 1 )(B E 2 )(B B)d 3 x, (47) 
where the coupling constant Ai = 13e 6 /3157r 2 m 8 and magnetic and electric parts of photon field 



are 

B = i{A>K) 1/2 k x E\\e- l{kr - ut) a\\{k), E x = -2(47r) 1/ V^±e i(feir -" lt) al(fc 1 ) (48) 

and similarly for E2, here uo = u)n(k) and ioi = u±(ki), i = 1,2. Here B, Bq, Ei are three- 
dimensional fields and ki are three dimensional vectors. 

For the decay || — > ||+ _L one has a similar interaction operator with the operator structure 

A+{t) = Aa+(p - k)a\_{k)a\\ (p)e~ itE (49) 

where E = oo»(p) — LUu(p — k) — oj±_(k). The coupling constant A in this case can be estimated 
as A/Ai = a(B /B cr ) 2 , where a is the fine structure constant, a = e 2 /he and B cr = m 2 c 3 /eh ~ 
4.4 x 10 13 Gauss. 

Let us consider a photon cascade created by a photon with momentum p and polarization ||. 
The photon splits as j\\(p) — > 7±(&i) + 7||(p — fci). Then one has the next generation of splitting: 
711 (p — k\) — > J±(k2) + 7|| (j> — k\ — k<i) etc. After generations of splitting one gets a cascade with 
N photons with _L polarization and momenta k\, k 2 ■ ■ ., k^ and also one photon with || polarization 
and momentum k — k\ — . . . — fcjv- An example of a cascade with two generations is shown in 
Fig.2. 

Our goal is to consider cascades with real photons (i.e. on the mass shell) including the 
intermediate states (compare with what one sees in the Wilson camera). These cascades can been 
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1 2 



2 1 



Figure 3: Non-crossing and crossing diagrams 

drawn as shown in Fig. 2. This diagram is not a Feynman one because all the lines (including an 
intermediate one) correspond to real particles on the mass shell and not to virtual states. More 
precisely all the lines in the diagram are "dressed" lines on the mass shell and the initial photon 
7||(fc) is prepared in a special way such that it undergoes the decay in a finite time. So we cannot 
use the standard S-matrix approach and the standard Feynman diagram technique to describe this 
process. The diagram is also not a diagram in the non-covariant diagram technique || because we 
have the conservation of energy at every vertex. The cascade in Fig. 2 may be intuitively described 
by the following state 

\iP{p, h, k 2 )) = f(p, h)f(p - h, k 2 )a+{p -h- A; 2 )a+(fc 2 )a+(A; 1 )|0) (50) 

where momentum conservation is built in creation and annihilation operators and energy con- 
servation is accounted for by the factor f(p, k) = f(u\\(p) — u>±(k) — u)\\(p — k)) where f{ui) is a 
function with support at u — 0. As we shall see below this is not the 5-function but roughly 
speaking its "square root". Indeed, the transition amplitude between two cascade states is given 
by scalar product 

(i/>(jpf, k[, k' 2 )m P , k h k 2 )) = \f(p, h)\ 2 \f(p - h, k 2 )\ 2 5(p - p'Wh - k[)5(k 2 - k' 2 ) (51) 

Notice that in the scalar product ( |5"T| ) only the non-crossing diagram Fig. 3a contributes. In 
fact the contribution from the crossing diagram in Fig. 3b is negligible because of conservation 
of energy and momentum. This is the crucial point where the difference between our diagrams 
describing real particles in intermediate states and the Feynman diagrams having virtual particles 
in intermediate states is manifested. In the Feynman diagram technique the amplitude of emission 
of the two photons is represented by a sum of two diagrams differing by the order in which the 
two photons are emitted. Here we have only one diagram, Fig. 2. 

Now let us observe that if in (|5H| ) we replace the operators a]_{ki) and a\_(k 2 ) by the quantum 
Boltzmann operators b±{ki) and bj_(k 2 ) satisfying the relations ([I]), i.e. b±(k)b±(p) = 5{k — p), 
then it will be automatically guaranteed that only the non-crossing diagrams survive. Therefore 
it is natural to describe cascade wave functions in terms of these operators. It is well known that 
standard free photons are bosons. Therefore to see the quantum Boltzmann statistics we have 
to prepare a special state depending on the interaction. In fact it is natural to expect that the 
cascades with physical intermediate states occur at a time scale slower than the one occurring in 
the standard S'-matrix approach to multiparticle production. A natural method, leading to this 
result, is suggested by the stochastic limit technique. 

Now let us consider the question how one can prepare a state which exhibits the new statistics 
for photons in cascade. If we would deal with the scattering of 2-particle states at infinite time (S'- 
matrix) we simply have to consider two Feynman diagrams to take into account the Bose statistics 
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of photons. However in the cascade we deal with evolution in finite time and the states of photons 
7_i_(fci) and j±(k 2 ) are prepared in a special way because they are emitted at times t\ and t 2 , 
respectively. Therefore, there is a reason not to add the second diagram. A special procedure 
which is adequate to this situation is the stochastic limit technique described in Sect. 3. In our 
case the master field is given by the asymptotic behaviour of the following collective operator 

lim -at(p - k)a+{k) ail {p)e- itE/x2 = B + {p, k) (52) 

where E is the same as in (f49|) and as in Sect. 3 - Sect. 5 the limit is meant in the sense of the 
Wightman correlation functions. 

As it follows from Theorem 1 the master field B (p, k) satisfies the following commutation 
relations 

B{p, k, t)B + (p', k', t') = 2nS(t - t')5{E)5{k - k')6(p - p')b+{k)b\\{k) (53) 
and 6||(p), 6j|"(p') satisfy the relations 

h(p) b t(p) = s (p~p') ( 54 ) 

The presence of the <5(£')-factor (E = uj\\{p) — oj\\{p — k) — u>±(k) ) has two important phys- 
ical consequences. First, the commutation relations for the B# are not a consequence of the 
corresponding relations for bj and b\: the three photons are entangled into a single new object 
(triphon). Second, the triphon creation and annihilation operators B# operate not on the usual 
Fock space but in interacting Fock space. 

By introducing the auxiliary creation and annihilation operators b(t),b + (t), b±(p),b±(p), sat- 
isfying the quantum Boltzmann relations 

b(t)b+(t')=5(t-t'), (55) 

b±(p)bl(p') = 5(p-p') (56) 
and introducing the symbolic relation 

B+(p,k,t) = b + (t)b+(k 1 )bi(k 2 )b ll (k)(27r) 1 / 2 5 1/2 (E) (57) 

we can disentangle the master field by expressing it as a product of individual Boltzmannian 
fields. Here the notation 6i/ 2 {E) is purely symbolic and it simply means that, since the new 
commutation relation (|53|) is quadratic in the master creation and annihilation operators, if in 
the above symbolic relation we consider the symbol Si/ 2 (E) as a scalar function satisfying the 
formal relation [5i/ 2 (E)] 2 = 6(E), and if we use the right hand side of (|57D to express the left 
hand side of ([53]), then the standard Boltzmannian relations (|56|) and (|54|) will reduce (|53|) to 
an identity. The intuitive understanding of the "disentangling" relation ( |57D is that the triphon 
master field B + (p,k,t) can be expressed as the product of three "Boltzmannian photons": one 
can think that each photon of the triple has its own (Boltzmannian) creation and annihilation 
operators depending on its own momentum, however in the master field these operators can only 
appear in the combination given by the right hand side of ( |57D and there is a constraint among 
the three momenta expressed by energy conservation. 

Notice the Boltzmannian white noise relation fl5"5|), which makes our model particularly suitable 
for Monte Carlo simulations. The origin for these new commutation relations lies in the fact that 
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the crossing diagrams in the computation of the matrix element ( pD| ) are suppressed in the weak 
coupling/large time limit. 

A photon splits into two not only in a magnetic field but also in a nonlinear medium. In 
fact such processes are well known in nonlinear quantum optics, see for example ||12|| . In the 



nonlinear process of parametric down conversion a high frequency photon splits into two photons 
with frequencies such that their sum equals that of the high-energy photon. The two photons 
produced in this process possess quantum correlations and have identical intensity fluctuations. 

In conclusion, in this paper we have argued that photon cascades in a strong magnetic field 
might create a new type of entangled states (triphons) which obey not Bose but the quantum 
Boltzmann statistics. This prediction is based on the assumption that both kinematically allowed 
photon splitting modes operate. 

Given the validity of this assumption we prove that, in the stochastic regime the interme- 
diate photons in a cascade are real and virtual particles. The dominating contributions to the 
dynamics come from these entangled triples of photons which behave like single new particles 
(triphons) whose statistics can be experimentally observed by counting the emitted photons in 
the corresponding cascades. 

As explained in the introduction, the time scale in which our predictions are true is long 
compared to the strength of the coupling but short if compared to the time scale of the S-matrix 
approach. This remark should be kept into account in a possible experimental verification of these 
predictions. A better theoretical understanding of the photon splitting with a non weak dispersion 
is required. From the experimental side new more precise devices such as the planned Integral 
mission [37] might significantly advance our understanding of the fundamental problem of photon 
statistics. 
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